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ABSTRACT 

We study the constraints which the next generation of radio telescopes could place 
on the nature of dark energy, dark matter and inflation by studying the gravitational 
lensing of high rcdshift 21 cm emission, and we compare with the constraints obtain- 
able from wide-angle surveys of galaxy lensing. If the rcionization epoch is effectively 
at z ~ 8 or later, very large amounts of cosmological information will be accessible 
to telescopes like SKA and LOFAR. We use simple characterizations of reionization 
history and of proposed telescope designs to investigate how well the two-dimensional 
convergence power spectrum, the three-dimensional matter power spectrum, the evo- 
lution of the linear growth function, and the standard cosmological parameters can be 
measured from radio data. The power spectra can be measured accurately over a wide 
range of wavenumbers at z ~ 2, and the evolution in the cosmic energy density can 
be probed from z ~ 0.5 to z ~ 7. This results in a characterization of the shape of the 
power spectra (i.e. of the nature of dark matter and of inflationary structure genera- 
tion) which is potentially more precise than that obtained from galaxy lensing surveys. 
On the other hand, the dark energy parameters in their conventional paramctrization 
(Aa, w o, w a ) are somewhat less well constrained by feasible 21 cm lensing surveys than 
by an all-sky galaxy lensing survey. This is because dark energy is felt primarily at 
relatively low redshifts in this model; 21 cm surveys would be more powerful than 
galaxy surveys for constraining models with "early" dark energy. Overall, the best 
constraints come from combining surveys of the two types. This results in extremely 
tight constraints on dark matter and inflation, and improves constraints on dark en- 
ergy, as judged by the standard figure of merit, by more than an order of magnitude 
over either survey alone. 

Key words: large-scale structure of Universe - dark matter - gravitational lensing 
- intergalactic medium - low frequency radio astronomy 



1 INTRODUCTION 



the energy density of the Universe is apparently dominated 
today by a different and even more unexpected field, dark 
energy, which is similarly accelerating the present cosmic ex- 
pansion. The nature of the dark energy and its relation to 
the rest of physics are unknown. Detailed measurements of 
the recent expansion history and of the corresponding gravi- 
tationally driven growth of structure provide the only known 
route to narrow down the possibilities. The techniques men- 
tioned above can measure the expansion precisely out to 
z ~ 4 and the growth of structure out to z ~ 1.5. In this 
paper we investigate a new technique that complements, re- 
fines and greatly increases the precision of these methods by 
directly constraining both cosmic expansion and structural 
growth out to z ~ 6. 



The now widely accepted cold dark matter (CDM) model 
combines with inflation-inspired scale-invariant primordial 
density fluctuations to provide a consistent explanation 
for cosmic microwave background (CMB) fluctuations, for 
type la supernova luminosity distances, for the clustering 
of galaxies in redshift surveys, for the galaxy cluster abun- 
dance and its evolution, and for the statistics both of weak 
gravitational lensing and of Lya forest absorption in quasar 
spectra. This success comes at the cost of introducing several 
mysterious and apparently ad hoc constituents. Most of the 
matter is supposed to be dark, a weakly interacting neutral 
particle so far detected only through its gravitational effects. 
All structure is supposed to have originated during a very 
early period of accelerated expansion driven by an inflaton 
field which has been posited purely for this purpose. Finally 



The spin temperature of neutral hydrogen during and 
before the epoch of reionization (8 ,$ z ,$ 300) fell out of 
thermal equilibrium with the CMB radiation, resulting in 



2 Metcalf & White 



the absorption and emission of 21 cm radiation. There has 
been a great deal of interest in the prospect of detecting and 
mapping this radiation using low-frequency radio telescopes, 
Several suitable instrum ents are now under con struction or 
in planning stages (see iFurlanetto et all 120061 . for an ex- 
tensive review). This radiation provides an excellent source 
for gravitational lensing studies. Structure is expected in 
the 21 cm emission down to arcsecond scales, and at each 
point on the sky there will be ~ 1000 statistically inde- 
pendent regions, separated in redshift (and thus frequency) 
that can in principle be observed. Gravitational lensing co- 
herently distorts the 21 cm brightness temperature maps at 
different frequencies. This coherent distortion can be dis- 
tinguished from intrinsic structure in the HI gas if enough 
independent redshifts are observed. In t his way a map of the 
foreground density can be constructed (|Zahn fc Z aldarriaga 
2009 ; iMetcalf fc Whitell2007l ; Iffilbert et alj|2007t iLu fc Pen 



20071 ). 



Observing the 21 cm radiation at high redshift will 
be challenging. It will require large arrays of radio tele- 
scopes working at low frequencies ( ~ 100 to 200 MHz). At 
these frequencies foregrounds - terrestrial, galactic and ex- 
tragalactic - are very large and will need to be subtracted by 
comp lex and as yet untested procedures (|Furlanetto et all 
120061 ). This challenge is currently being addressed by sev- 
eral observational teams (see section [6j|. 

In addition to the noise associated with mapping the 
brightness temperature, the lensing signal has an additional 
intrinsic noise component which comes from the unknown 
intrinsic structure of the 21 cm brightness temperature 
distribution. This noise cannot be reduced by increasing 
the collecting area of the telescope, by increasing the in- 
tegration time or by im proving the removal of foregrounds. 
IMetcalf fc White! (|200?n showed that if the signal-to-noise in 
the brightness temperature map at each frequency is greater 
than one, then the noise in the mass map will be close to 
the intrinsic value. Increasing the frequency resolution of the 
radio observations increases the number of effectively inde- 
pendent regions along the line-of-sight until the bandwidth 
becomes smaller than the correlation length in the redshift 
direction of the brightness temperature distribution. If the 
bandwidth is matched to the correlation length, the intrinsic 
noise is minimized. The correlation length in turn depends 
on beam size, and is smaller for smaller beams. Thus unlike 
galaxy lensing surveys, the intrinsic noise for 21 cm lensing 
decreases with increasing telescope resolution. In practice, 
there is a trade-off because smaller bandwidth means less 
flux, but this can be compensated by increasing collecting 
area and/or integration time. The noise is also affected by 
the range of frequency (i.e. redshift) over which the 21 cm 
radiation can be detected. The low-frequency limit is set 
by the telescope and/or the ability to subtract foregrounds. 
The high-frequency limit is typically set by the reionization 
of the universe, after which the amount of intergalactic HI 
is small. 

In lHilbert et~al] (|2007l ) we simulated how well an ideal- 
ized optimal telescope would be able to map the 2 dimen- 
sional distribution of matter. In this paper, we study how 
well radio telescopes with specifications similar to those of 
instruments currently being built or planned will be able to 
constrain cosmological parameters, in particular those re- 



lated to the nature of dark matter, dark energy and infla- 
tion. 

This paper is organized as follows. The next two sec- 
tions introduce the formalism we use to study lensing in 
general and lensing of the 21 cm radiation in particular. In 
section [3] we develop a formalism for extracting cosmologi- 
cal information and density maps from such lensing data. In 
section [5] we discuss the relevant aspects of reionization and 
the simple reionization model we use to make quantitative 
predictions. The telescope parameters used in our predic- 
tions are given in section [6] Section [7] contains predictions 
for the noise levels in various quantities. A summary of our 
results and of future prospects is given in the last section. 



2 LENSING PRELIMINARIES 

Gravitational lensing shifts the observed position of each 
point in the image of a distant source. Take the observed 
angular position on the sky to be 6 and the position in the 
absence of lensing to be j3. The first-order distortion in the 
image is expressed by the derivatives of the mapping be- 
tween these angles. The distortion matrix is commonly de- 
composed into the convergence k and two components of 
shear, 7, defined by 
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To lowest orde r and to an excellent approximation 
{Vale fc White! 120031 ) the convergence is related directly to 
the distribution of matter through 
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The weighting function for the source distance distribution, 
r](z), is normalized to unity. D(z' , z) is the angular size dis- 
tance between the two redshifts and 5(x, z) is the fractional 
density fluctuation at redshift z and perpendicular position 
x. The function 



e{z) = yUn(i + zY + n A (i + zyuz) + (] 



■Oa)(1-t 
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where fl m and Qa, are the present day densities of matter 
and dark energy measured in units of the critical density. 
The function describing the evolution of dark energy with 
redshift can be written 
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where w(a) is the equation of state parameter for the dark 
energy - the ratio of the of its pressure to its density - and 
a=(l + 2) _1 is the scale parameter. Where not otherwise 
mentioned, we will assume the universe is flat, Q m +Qa = 1, 
and w(a) — — 1. 
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Equation Q shows that k can be interpreted as a kind 
of projected dimensionless surface density. For our purposes 
it is convenient to express equation ([3]) as a matrix equation, 



K = G<5, 



(7) 



where the the components of K are the convergences run- 
ning over all position angles 9 and source redshifts, z s . The 
components of the vector 8 run over all position angles and 
foreground redshifts Zi. The matrix G is a function of most 
of the global cosmological parameters - Q m , ^a, w, etc - and 
is independent of position on the sky. The latter property 
makes these equations equally valid when k and 8 are trans- 
formed from angular space to spherical harmonic space or to 
the u-v plane where interferometer observations are carried 
out. 

When considering 21 cm lensing we will make the ap- 
proximation r](z, z s ) = 8 D (z — z s ), a Dirac delta function, 
for each band observed corresponding to a frequency of 
v = 1420(1 + Zs)^ 1 MHz. This is reasonable because an- 
gular size distances vary little within a single band, but it 
is not true in the case of galaxy lensing which we address in 
section 17.41 



3 LENSING OF PREGALACTIC HI 

Many convergence estimators are possible and can be ex- 
pressed either in real-space (which is more easily visualised) 
or in Fourier-space (which is more easily related to inter- 
ferom eter observations in visibility-space ). iMetcalf fe Whitel 
120071 ) gave a real-space estimator for k that has some ad- 
vantages, but to simplify the present exposition we here 
use a Fourier-space es timator based on that introduced by 
iHu fc Okamotol (|2002|) for the case of CMB lensing, and ex 



te nded by Zahn fc Zaldarriagal |2006h (and in the appendix 
of IMetcalf fc Whitel (120071 )1 to allow estimation of the 2- 
dimensional k field from 3-dimensional 21 cm data. Our 
treatment here parallels this earlier work, although with 
some significant differe nces in the source redshi ft weighting. 

In the appendix of IMetcalf fc Whitel (|2007n ) we applied 
the IHu fc Okamotol (|2002h estimator for the Fourier trans- 
form of the convergence, k(£, u), to each frequency band and 
then combined them optimally to find a final estimate for 
k(£). This has the disadvantage of requiring that an optimal 
bandwidth be found numerically which will be a function of 
both \£\ and frequency. An alt ernative method for combin - 
ing frequencies is developed in IZahn fc Zaldarriagal [|2006h . 
They use the Fourier transform of T(£, v) in the frequency 
direction as well as the angular directions and then find a 
convergence estimator of the form 

k(l)= f d 2 t'Y^Y. %{«■', t-MM)^ {£', k^)T* {£' -£,k 2 ) 

J fc x k 2 

(8) 

where £ is the Fourier dual of the angle and k is the 
discrete Fourier mode in the radial direction. The kernel 
x(£',£, ki, fe) is found by minimizing the noise while con- 
straining the average to be the convergence. This has the 
apparent advantage that the bandwidth does not come into 
the calculation, but, in practice, an equally arbitrary fre- 
quency scale will need to be imposed. In deriving the noise 



in their estimator IZahn fc Zaldarriagal (|2006l ) make the im- 
plicit assumption that the temperature field is statistically 
homogeneous in frequency. This assumption results in the 
Fourier modes being uncorrelated which greatly simplifies 
the calculation. In practice the temperature field will be in- 
homogeneous because the noise will be a strong function of 
frequency and because of the evolution of structure and ion- 
ization. As a result their estimator will not be optimal for a 
finite range in frequency only for Ai//u <C 1. This problem 
was avoided in our method. 

In this paper we will adopt a hybrid method for estimat- 
ing the convergence that avoids many of the shortcomings 
the pr evious two methods. We use the lZahn fc Zaldarriagal 
( (20061 ) estimator within frequency bands which are larger 
than our previous bandwidths and the minimum band width 
of the telescope, but small enough that no significant inho- 
mogeneity in the noise and brightness temperature structure 
is expected. This avoids having to find the optimal band- 
width. We then linearly combine these n(£, v) estimates in 
optimal ways as described in detail in section [4] We show 
that the correlations between convergence estimates at dif- 
ferent frequencies is not significant if the bands widths are 
taken to be ~ 1 MHz. 

The optimal kernel for estimator JS} under the assump- 
tion of homogeneity and Gaussianity within the band is 

X(£',£,ki,k2,v) - u{t, v)\£\ 2 8 klk2 

[£ ■ £' C„{£' M) +£■ {£ -£')C V {\£' - £\M)\ 
X C?(£>,ki)CT(\£' -£\,ki) 

(9) 

where C„(£, k) is the power spectrum of the actual temper- 
ature, while C v {Z,k) = Cj(£, k) + C„{£,k) is the observed 
power spectrum which includes noise. The normalization is 
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In the limit of an infinitely large beam the correlation 
between modes is 

(k (L, fci) k* (L', fe)> = 2(2ty) 4 8 d (L - L') 8 klk2 

x J d 2 £' x(^',L,fei) 2 Cj(/,fei)Cj(|£' -L|,Jfei) (11) 

= (2n) 2 S D (L - L') S klk2 N k (L, v) (12) 

To simplify analysis, it has here been assumed that the tem- 
perature (i.e. the 21 cm emissivity) is Gaussian distributed 
so that the fourth moment can be written as products of 
second moments. The actual temperature distribution will 
be non-Gaussian because of non-uniform spin temperature, 
peculiar velocities, nonlinear structure formation, and non- 
uniform ionization. The importance of nonlinea r structure 
forma tion in this context has been highlighted bv lLu fc Pen! 
(2007). Some of these non-Gaussian effects, such as non- 
uniform ionization, cannot be usefully quantified at this 
time. The Gaussian approximation should be valid at least 
before significant reionization occurs and is probably good 
while the reionized regions are significantly smaller than the 
resolution of the telescope. 

The finite telescope beam will cause correlations in the 
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noise, {k(£)k* (£ + Si)) 7^ 0, when 5 £ ^ 2tt<j u where a u is th e 
width of the beam in u-v space (|MetcaIf fc White! [2007h . 
These correlations can be taken into account, but for sim- 
plicity they will not be considered here. 

The noise in k{£, v) within one frequency band is 



N K (£,u) 

?/ 



(2^) 2 



d 2 £' 



[£ ■ £' C v (£\ k) +£■(£- 0) C v (\i' - £\, k)f 



C?(£',k)CT(\£' -£\,k) 



(13) 

Because the estimator JH} is a sum over all the observed 
pairs of visibilities, it will (by the central limit theorem) be 
close to Gaussian distributed even though it is quadratic in 
the visibilities. In the remainder of this paper this property 
will be assumed. The intrinsic noise limit corresponds to the 
case where C v {£) = C%(£). 

In real data the temperature distribution will not be 
Gaussian, the foregrounds will not be perfectly subtracted 
and will produce spurious frequency correlations, there will 
be holes in the area surveyed around bright point sources, 
there will be a finite and irregular beam, and the coverage 
of the u-v plane will not be complete. It is also true that the 
convergence estimator was derived in the weak lensing limit 
(k <C 1) which will not be valid in all regions of the sky. All 
these complications make it unclear at the present time what 
estimator will be optimal for a real experiment. We never- 
theless believe that the above relatively simple assumptions 
should give a good indication of what can be expected. 



4 PARAMETER ESTIMATION AND DENSITY 
MAPPING 

The noise in k(£, z) for different bands and different I will, 
to a good approximation, be statistically independent. The 
data vector will be defined as 



D = K K 

= K G5 



(14) 
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where the components run over all the combinations of Zi 
and £ that are measured. 

In this case the log of the likelihood function can be 
written 



ln£ 
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An additional function H has been added to represent a 
prior distribution on the parameters or, in the context of 
density reconstruction, a regularization. As written so far 
the free parameters include all the cosmological parameters 
and all the foreground densities in each redshift bin. The 
noise covariance matrix is approximately 



^5 ij I(£ i ,v i )N ii (£i,n). 



(17) 
(18) 



The matrix J, y „/(L) has been introduced to express pos- 
sible cross-correlations between frequency bands. It will 
be normalized so that 7„„(L) = 1. The second line fur- 
ther simplifies this by assuming that the noise is equal in 
frequency bands that have significant correlation between 



them. In is case the I vv i (L) matrix can be made into a fac- 
tor, I(£, v) = Y^ v i Ivv' (£)■ In our studies here I vv i (L) ~ 5 vv i 
or I(£, v) ~ 1 since the correlations between the wide bands 
used (~ 1 MHz) are small, but it is possible that fore- 
ground subtraction in particular might introduce significant 
correlations between frequency bins. Foreground subtraction 
will not be discussed in detail here, but we will retain the 
off-diagonal elements in our formalism. There will also be 
some off-diagonal elements to Nij between different £ values 
caused by the finite beam of the telescope. Thses could be 
incorporated in a future analysis. 

The maximum likelihood estimate for any parameter 
can be found by maximizing (|16[) with respect to that pa- 
rameter. The error in this estimator is often forecasted using 
the Fisher matrix defined as 



d 2 \nC 

dpidpj 



(19) 



The expected error in the parameter p a , marginalized over 
all other parameters, is a 2 ~ (F" 1 )^. The unmarginalized 
error estimate (the error when all other parameters are held 
fixed) is (F aa y\ 

Often the parameters of physical interest have highly 
correlated noise, quantified by the off-diagonal elements in 
F. As a result a a can be a deceptive measure of how well 
the data constrains the parameter set as a whole. A way to 
mitigate this is to find the transformation that diagonalizes 
F 



V f AV. 



(20) 



This defines linear combinations of the parameters p = Vp 
that are uncorrelated and have variances 

Aoo ■ This will be 

used in section [7731 



4.1 Tomography 

If one is primarily interested in reconstructing the cosmic 
mass density distribution, the background cosmology can 
be held fixed and C can be maximized with respect to the 
pixelized foreground density. This can be done in 2-D by 
approximating re as independent of v, in which case the so- 
lution that maximizes (|16|1 is 
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N(£) 
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(22) 
(23) 



which is uncorrelated between £ values that are separated 
by more than the resolution of the telescope. Here no regu- 
lar izing function is used (H — 0). This is the result derived 
bv lMetcalf fc White! (|2007h . 

A 3-D density reconstruction can be obtained by maxi- 
mizing (|16[) with respect to the components of S. For H — 0, 
the solution, after some algebra and taking into account that 
N is symmetric, is 



g^^gI G+N^K. 



(24) 
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Essen tially the same result was obtained by iHu fe Keetonl 
(12002) 1 for density reconstruction using galaxy lensing data 
in real-space. This requires that the matrix in braces be 
invertible. In the special case where G is invertible the es- 
timator is simply <5 = G _1 K. The noise covariance matrix 
for the reconstruction ()24p is 



N 4 



[G + I 



(25) 



This noise will be highly correlated between 8's. 

The regularization function in (|16[1 can be used to im- 
prove the density reconstruction at the expense of making 
some assumptions about the statistical properties of the un- 
derlying density distribution. One choice is to assume that 
this distribution is Gaussian, in which case 



5{z t ,£f 
Cs(zi,i)' 



(26) 
(27) 



where Cs{z, €) is the angular power spectrum of mass within 
the redshift bin labeled by Zi. This is equivalent to a Wiener 
filter. The optimal estimator Q24| l an d the noise ([25} will 
be modified in this case, but can be derived in the same 
way. In addition to smoothing the noisy map, regulariza- 
tion provides a well behaved way to deal with missing tele- 
scope baselines and with holes in coverage due to foreground 
sources. Other regularization schemes could include an en- 
tropic prior or filters desig ned to emphasize localized mass 
lumps <|Hu fe Keetonll2002l ). 



4.2 Cosmological parameter estimation 

Estimates of cosmological parameters and statistical infor- 
mation about the underlying mass density distribution can 
be extracted from 21 cm data without constructing density 
maps. The actual distribution of matter can be marginalized 
over assuming a suitable prior. This procedure can give sur- 
prisingly precise results, even when noise levels are far too 
high for a meaningful density reconstruction to be possible. 

According to the standard model of structure forma- 
tion, the components of K (i.e. v)) w ill be normally 
distri buted for £ less than several thousand (jTakada fe Jainl 
1 2004| ). In this case a prior H = ^K t C^ 1 K can be used 
and the likelihood function can be integrated over all com- 
ponents of K - i.e. over all possible convergence maps. C K 
here is the (cross-)power spectrum of the convergence for 
two different source redshifts. 



(28) 



This can be calculated using expression ([2]) and a model for 
the matter power spectrum. The resulting likelihood func- 
tion has the same form as (|16|) 



ln£ = -iK t C- 1 K-i|C|, 



(29) 



where the noise matrix, N, is now replaced with the covari- 
ance matrix, 

C = N + C K . (30) 
Since £-modes separated by more than the resolution of 



the telescope will not be correlated, we can break the likeli- 
hood function u p into factors representi ng each resolved re- 
gion in ^-space (|Metcalf fe White! [20071 ). The result is that 
there are ~ (2£ + l)/sky independent measured modes for 
each value of £, where / s ky is the fraction of sky surveyed. 
The Fisher matrix can then be further simplified to the 
widely used form, 



F ab = \ {21 + l)/ sky tr [C^C,. C^Gs, 



(31) 



A quantity of particular interest that will be calculated 
later in this paper is the noise variance in the 2-D conver- 
gence power spectrum estimate. This is easily derived from 
(f3"T|) as 



AC K {£) = 



(2£ + l)f aky 



\C K {t)+N{l)), 



(32) 



where N{£) is given by (|22[) . The first term represents the 
sample or cosmic variance, and the second the noise in the 
k estimate itself. To increase the signal-to-noise the average 
power (known as the band-power) can be estimated within 
a bin of width A£ larger than the resolution limit of the 
survey, S£ ~ fs\^/ 2 - The noise in the band-power is the above 
divided by the square root of the number of independent 
measurements within the band, 



AC K (£, A£) 



ACJfi 



(33) 



5 REIONIZATION MODEL 

The fluctuations in the brightness temperature depend on 
the spin temperature, T s , the ionization fraction, ih and the 
density of HI through 



5T b ~ 24(1 + S b )xn 
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(|Fieldl Il959l ; iMadau et all 119971 ). As is commonly done, we 
will assume that the spin temperature is much greater than 
the CMB temperature. This leaves fluctuations in xh, and 
the baryon density S b — (p b —~p b )/~p b as the sources of bright- 
ness fluctuations. We will make the simplifying assumption 
that xh = 1 until the universe is very rapidly and uniformly 
reionized at a redshift of z TC i on . Realistically, the reionization 
process will be inhomogeneous and may extend over a sig- 
nificant redshift range. This will increase C v (€) by perhaps 
a factor of 10 o n scales larger the chara cteristic size of the 
ionized bubbles (jZaldarriaga et al"1l2004l ) and thus might be 
expected to reduce the noise in the lensing map, k(9), signif- 
icantly. However, the noise, (|13p . has been derived assuming 
that the fourth order statistics of 8T b can be approximated 
by the values appropriate for a Gaussian random field. If 
this approximation holds, the lensing noise will indeed be 
reduced during reionization, but this is uncertain, since the 
temperature distribution will clearly not be Gaussian at this 
time, especially when the neutral fraction is low. Incorporat- 
ing a more realistic ionization history is difficult and requires 
further work. A definitive resolution of these uncertainties 
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will not, in any case, be possible until the observations are 
taken. 

With real data one will not need to rely on an assump- 
tion of Gaussianity. It will be possible to derive the un- 
certainties by bootstrap techniques, since the higher order 
statistics of the temperature variations can be estimated di- 
rectly from the data themselves. 



6 MODEL TELESCOPES 

In order to forecast the capabilities of future instruments, 
we will present results assuming several sets of telescope pa- 
rameters. These are intended to represent the characteristics 
of facilities currently under construction or in an advanced 
stage of planning. 

The noise in each visibility measurement will have a 
thermal component and a component resulting from imper- 
fect foreground subtraction. Here we model only the ther- 
mal component. If the telescopes in the array are uniformly 
distributed on the ground, the average integration time for 
each baseline will be the same and the power spectrum of 
the noise will be 



r*™ — 

w — 



2tt 



(2tt) 3 T 2 



/ cover max 



(35) 



(IZaldarriaga et aL | 12004 iMoralesI 120051 : iMcQuinn et~al] 
l2006h where T„, a is the system temperature, Av is the band- 
width, t is the total observation time, D te i is the diameter 
of the array and ^ max (A) = 2nD te \/\ is the highest multi- 
pole that can be measured by the array, as set by the largest 
baselines. / CO ver is the total collecting area of the telescopes 
divided by 7r(D te i/2) 2 , the aperture covering fraction. Other 
telescope configurations are possible which would result in 
the noise being unequally distributed in £, but here we will 
consider only this uniform configuration. At the relevant fre- 
quencies, the overall system temperature is expected to be 
dominated by galactic synchrotron radiation. We will ap- 
proximate the brightness temperature of this foreground as 
T sky = 180 K(z//180 MHz)" 2 ' 6 , as appropriate for region s 
well away from the Galactic Plane (|Furlanetto et al1l200rj) . 
This results in larger effective noise for higher redshift mea- 
surements of the 21 cm emission. 

Several relevant telescopes are currently proposed or 
under construction. The 21 Centimeter Array (21CMA, for- 
merly known as PAST)Q has /cover ~ 0.01 and ^ max ~ 10 
giving it a resolution of about 10 arcmin. The Mileura Wide- 
field Array (MWA) Low Frequency Demonstrator (LFDj3 
will operate in the 80-300 MHz range with D tc i — 1.5 km 
and / covcr ~ 0.005. For LOFAR (the Low Frequency Array jf| 
the core array was originally planned to have / covcr ~ 0.016 
and Dtc\ ~ 2 km. LOFAR's extended baselines, reaching 
to 350 km and possibly farther, will not be useful for high 
redshift 21 cm observations because of the small / COVGr , al- 
though they will be helpful for foreground subtraction. It is 
anticipated that LOFAR will be able to detect 21 cm emis- 
sion out to redshift z ~ 11.5, but sensitivity limitations will 
make mapping very difficult. Recently budget constraints 



1 21cma.bao.ac.cn/ 

2 www.haystack.mit.edu/ast/arrays/mwa/ 

3 www.lofar.org 



have reduced the initial collecting of the telescope core by 
about one third. We will use approximations to these orig- 
inal plans for our forecasting and call our mock telescope 
LOFARII to indicate that it is larger than the first stage of 
LOFAR being constructed. Plans for the SKA (the Square 
Kilometer Array )Q have not been finalized, but it is expected 
to have /cover ~ 0.02 over a diameter of ~ 6 km (£ ma x ~ 10 ) 
and sparse coverage extending out to 1,000-3,000 km. The 
lowest frequency currently anticipated is ~ 100 MHz which 
corresponds to z ~ 13. It is expected that the core will 
be able to map the 21 cm emission with a resolution of 
A8 ~ 1 arcmin. For reference, one arcminute (fwhm) corre- 
sponds to baselines of 5.8 km at z — 7 and 11 km at z = 15. 

The various sets of telescope and survey parameters we 
use below are listed in Table [T] We adopt values appropriate 
to one or three seasons of observing with the core arrays of 
LOFARII and SKA, and we assume that reionization occurs 
instantaneously at z = 7 or (for one SKA case) at z = 10. 



7 FORECASTS 

In the following subsections the formalism of sections [5] 
through [4] and the models of sections [5] and [6] are used 
to forecast how well future projects will be able to esti- 
mate various quantities of interest. As background cosmo- 
logical model we choose f2 ma tter = 0.3 (matter density), 
= 0.7 (cosmological constant), H = 70 (Hubble param- 
eter), as = 0.75 (normalization of the matter power spec- 
trum) and Q,b = 0.03 (baryon density). Except in section [731 
the (nonlinear ) matter power spectrum is calculated using 
the method of lPeacock fc Doddsl fl996). In the following we 
will make the approximation I(L, v) = 1 which will be fur- 
ther justified. 



7.1 Convergence power spectrum 

The first objective of a 21 cm lensing survey might be to 
measure the two-point statistics of the convergence field 
k(£, z s ), averaged over z 3 . The expected error in the binned 
power spectrum C K is given by equations (|32[) and ((33} . For 
purely aesthetic reasons we choose the ^-space bin-widths to 
be 



Al = 



/bin 

/sky 



C K (^=10) AC K (£) 



AC K (£ = W) C K (£) 



10 



3/5 



(36) 



Here, the bins at £ = 10 are /bin times the telescope resolu- 
tion, where the factors /bin used for figures [T] through [4] are 
given in the captions. 

Figure [1] shows forecasts for a LOFARII- like telescope 
with observing times of 30 days (approximately one season) 
and 90 days. Reionization is assumed to occur at z = 7. 
It can be seen that the noise per mode is nowhere smaller 
than the signal per mode, even after three seasons. Thus, 
LOFARII will be unable to make a high-fidelity map of the 
density distribution unless reionization is very different than 
assumed here. Despite this, after three seasons the band- 
averaged convergence power spectrum can be recovered with 
good signal-to-noise over a wide range of £- values. 



www . skatelescope . org/ 



Cosmological Information in the Gravitational Lensing of Pregalactic HI 7 



Table 1. Parameters for telescope models. 



model 


D (km) 


f 'cover 


t (days) 


!/ min (MHz) 


Av (MHz) 


^rcion 


rcionization 


LOFARII-lyr 


2 


0.016 


30 


110 


1 


7 


instantaneous 


LOFARII-3yr 


2 


0.016 


90 


110 


1 


7 


instantaneous 


SKA-lyr 


6 


0.020 


30 


100 


1 


7 


instantaneous 


SKA-3yr 


6 


0.020 


90 


100 


1 


7 


instantaneous 


SKA-3yrA 


6 


0.020 


90 


100 


1 


10 


instantaneous 





Figure 1. Forecasts of the lcr uncertainties in estimates of the 

convergence power spectrum for our LOFARII-lyr (blue) and 

LOFARII-3yr (red) parameter sets. In the upper panel the solid 

straight lines give the noise per mode, while the dotted curve 

is the underlying model power spectrum. We have here assumed 

— 1/2 

that 25% of the sky has been observed (the noise scales as / aky )• 
The band powers are for bins in i chosen according to formula J36H 
with /bi n = 10. The signal-to-noise ratio for the banned power is 
given in the lower panel. 



Figure [2] shows similar forecasts for an SKA-like tele- 
scope. Here the noise per mode is below the signal for 
I < 700 after 30 days observing, and for i < 2000 af- 
ter 90 days, so mapping the convergence on these scales 
should be possib l e. Th is agrees with the conclusions in 
iMetcalf fc White! (|2007t) . although we used a different cri- 
terion for judging a "good" map in that paper. In cases 
such as this, where the noise per mode is smaller than the 
expected power in the signal, the uncertainty in an estimate 
of the power spectrum is dominated by cosmic variance, and 
the only way to increase the precision of the measurement 
is to survey a larger fraction of the sky. This can be seen 
in figure [5] where the noise becomes independent of observ- 
ing time at small i. With the parameters assumed here, an 
SKA-like telescope will clearly be able to make an excellent 
estimate of the convergence power spectrum. 

Figure [3] illustrates how uncertainties in the power spec- 
trum estimate depend on the redshift of reionization. We 
compare results for instantaneous reionization at redshifts 
of 7 and 10 for 3 years of observation with an SKA-like 
telescope. When z rc i on is high, there are fewer independent 
sources of 21 cm emission in the observable frequency range 
and the noise goes up significantly. Nevertheless, using only 
the emission between z = 10 and 13, a good convergence 
power spectrum can still be recovered, although mapping 



Figure 2. Forecasts of the lcr uncertainties in estimates of the 
convergence power spectrum for our SKA-lyr (blue) and SKA-3yr 
(red) parameter sets. The band powers are for bins in I chosen 
according to formula (136 I t with /t,i n = 5. All other parameters 
are the same as in Figure [T] Note that below about I ~ 10 3 
the uncertainties are dominated by cosmic variance and do not 
improve for the longer observing time. 




Figure 3. Forecasts of the lcr uncertainties in estimates of the 
convergence power spectrum for our SKA-3yr (red) and SKA- 
3yrA (blue) parameter sets. The band powers are for bins in I 
chosen according to formula (136 I t with /t,; n = 10. All other pa- 
rameters are the same as in Figure[l] The difference between these 
two examples is the reionization redshift, 2reicm — 

7 for SKA-3yr 

and z r eion = 10 for SKA-3yrA. 
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Figure 4. An illustration that increasing the frequency band- 
width of the observations results in increased noise in the k esti- 
mate. Red shows results for our SKA-3yr parameters which adopt 
a bandwidth of 1 MHz, while blue adopts the same parameters 
except that the bandwidth is taken to be 5 MHz. The band pow- 
ers here are for t bins with /bin = 10. All other parameters are 
the same as in Figure [T] 

the dark matter distribution would no longer be possible. A 
telescope that goes to lower frequencies could tolerate higher 
reionization redshifts without such performance degrada- 
tion. An example is the MWA, although it currently has 
insufficient collecting area for this project. A LOFAR-like 
telescope would be unable to obtain a useful power spec- 
trum estimate for z rc ion ^ 9. 

Figure [4] illustrates how the uncertainties depend on the 
bandwidth. Increasing the bandwidth by a factor of 5 makes 
little difference to the results. It is assumed here that the 
statistical properties of the noise and source 21 cm radiation 
are constant within each band. This result indicates that the 
correlations between frequency bands when 8v = 1 MHz are 
not of relevance since they are fully incorporated when the 
bandwidth is increased. As discussed before, the estimator is 
optimized for the case where the statistical properties of the 
noise and source 21 cm radiation are constant across each 
band and so it is better to use a smaller bandwidth where 
this assumption is better justified. 

The lensing map of the cosmic mass distribution can 
have good fidelity while the temperature map of high red- 
shift 21cm emission is noise-dominated on the same scale. 
This somewhat counter-intuitive situation reflects the fact 
that it is better to have more independent redshift slices at 
low signal-to-noise than to make high-quality images of the 
brightness temperature in a small number of channels. 

7.2 Matter power spectrum 

Two-dimensional k(tj measurements can also be used to es- 
timate the power spectrum of matter fluctuations directly, 
if the background cosmology (and thus the evolution with 
redshift of the matter fluctuations) is assumed to be known. 
To show this, we model the matter power spectrum by lin- 
early interpolating between points that are evenly spaced in 
In k. The In P(k) values are treated as the model parameters 
in this case. Equation (|31[) can then be used to find the ex- 




0.001 0.010 0.100 1.000 10.000 100.000 

k 

Figure 5. Forecast uncertainties in a reconstruction of the matter 
power spectrum from lensing of high-redshift 21 cm emission for 
observations with our LOFARII-3yr parameter set. We have taken 
7.93 points per decade in the model power spectrum and have 
used no tomographic information. We assume f s k y = 0.25. 




0.001 0.010 0.100 1.000 10.000 100.000 

k 

Figure 6. Forecast uncertainties in a reconstruction of the matter 
power spectrum from lensing of high-redshift 21 cm emission for 
observations with our SKA-3yr parameter set. We have taken 
15.87 points per decade in the model power spectrum and have 
used no tomographic information. We assume / s i; y = 0.25. 

pected uncertainties in the power spectrum reconstruction. 
For this latter calculation we assume purely linear evolu- 
tion with redshift and fix the shape of the pow er spectrum 
to be that given by the CAMB computer code (|Lewis et al.l 

l2oooh Fl 

The uncertainties in these P(k) measurements at differ- 
ent k values will be correlated, and their values will depend 
on the number of points used to represent the power spec- 
trum - the more points the better the fc-space resolution, but 
the larger the noise per point. Figures [S] and [5] show forecast 
uncertainties for two of our observational parameter sets. 
The matter power spectrum should be extremely well de- 
termined with 3 years of SKA data, and using tomographic 
information would improve the estimate still further. 

5 http:/camb.info 
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Figure 7. The uncorrelated modes in a decomposition of an 
estimate of the structure growth function for the LOFARII-3yr 
case. The upper panel shows the redshift dependence of the two 
modes which are expected to yield estimates with signal-to-noise 
above 1, while the lower panel gives the signal-to-noise values 
forecast for the 10 best constrained modes. The x-axis is simply 
the mode number. Corresponding modes are indicated using the 
same colour in each panel. The growth function is discretized into 
20 points spaced evenly in log(^) and k(£, z s ) are binned into 20 
evenly spaced z s bins. We assume / sky = 0.25. 

The comoving scales which contribute most to 21 cm 
lensing can be seen clearly in figures [5] and [6] This technique 
is complementary to CMB fluctuation measur ements which 
prob e structure on scales with 10 -4 < k < 0.1 (|Spergel et alj 
120071 ). The scales probed by 21 cm lensing and by galaxy 
clustering surveys overlap partly, but in the former case the 
mass is probed directly with no bias uncertainties, and the 
effective redshift of the measurement is much higher (see the 
next section) so that uncertainties due to nonlinearities are 
also reduced. Baryon acoustic oscillations can be seen in the 
power spectrum in figures [5] and [6] in the range k — 0.04 - 
0.4 Mpc -1 . They would be well measured in the SKA-3yr 
case. These features would provide additional leverage when 
estimating cosmological parameters, but we do not take this 
into account in section \7A\ below. 

7.3 Tomographic information 

As discussed in section |4j there is tomographic information 
in the measured k(£, z s ). This makes it possible, in principle, 
to measure how P(k) evolves with redshift, and so to probe 
dark energy through its indirect effects on the growth rate 
of linear fluctuations. 

While density perturbations are small they grow lin- 
early, 8(1, z) = D g (z)S(£, z = Zo) I 'Dg(zo) where z is some 



Figure 8. The uncorrelated modes in a decomposition of an es- 
timate of the structure growth function for the SKA-3yr case. 
The upper panel shows the redshift dependence of the five modes 
which are expected to yield estimates with signal-to-noise above 
1, while the lower panel gives the signal-to-noise values forecast 
for the 10 best constrained modes. Corresponding modes are indi- 
cated using the same colour in each panel. The growth function is 
discretized into 50 points spaced evenly in log(^) and the k(£, z s ) 
are binned into 20 evenly spaced z s bins. We assume f s ^ y = 0.25. 

initial redshift. Within General Relativity the linear growth 
function D g (z) is directly related to the cosmic expansion 
history a(z) and so through the Friedmann equation to the 
history of the various contributions to the cosmic energy 
density. In other theories of gravity D g (z) a nd a(z) are in- 
dependent (jPeeblesll 1980l ; lBertschingerll2006l ) . A comparison 
of precise measurements of the two can thus be used to test 
Einstein's theory. 

Although structure growth is not precisely linear for all 
the i values and for the full range of redshifts probed by 
21 cm lensing, it is probably a good approximation for the 
modes that are probed with high signal to noise. Without a 
specific theory for the nature and amount of dark energy the 
exact form of D g (z) is unknown. In this section we model 
D g (z) as piecewise linear between a set of interpolation red- 
shifts at which we treat its values as unknown parameters to 
be estimated. Note that this implies more freedom than is 
physically allowed, since D g (z) should, for example, decrease 
smoothly and monotonically with redshift. Technically our 
approach is not fully self-consistent since the angular size 
distances in k(£, z s ) [equation p|] depend on the cosmic ex- 
pansion history and so should change when D g (z) changes. 
Such effects are weaker than those coming from the growth 
function itself, so we neglect them here. In section PTTil below. 
we investigate an alternative approach in which we param- 
eterize the equation of state of dark energy and then use 
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Table 2. Forecast constraints on cosmological parameters for full sky surveys. All constraints scale with sky coverage as / aky except 
that the dark energy figure of merit, FOM, scales as / s ]jy. The quantities in parentheses are the marginalized uncertainties for constant 
w (w a = 0). Both the high redshift 21 cm emission and the lensed galaxy population are assumed to have been binned into ten disjoint 
redshift intervals when making the corresponding parameter estimates. 







In A 


n s 


r 


n A 


w 


w a 


FOM 




assumed values 


-y.D [erg — U. to) 
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LOFARII-3yr 


marginalized 
unmarginalizcd 


2.2 
n no/i 7 


0.22 
U.UUzl 


0.11 
n nnnQC 

U.UUUoo 


0.64 

U.UUZO 


3.0 (1.3) 

U.UZ ( 


11 
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0.075 


LOFARII-3yr 
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0.036 


0.0036 


0.0012 


0.0060 


0.43 (0.053) 


1.6 


12 


1 TT> 1 o n rAr 
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U.UU4y 


U.UU1D 
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U.UUOO 
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LOFARII-3yr 


marginalized 


0.032 


0.0045 


0.0017 


0.0033 


0.028 (0.0087) 


0.087 


1300 


+ galaxies 


unmarginalizcd 


0.00032 


0.00030 


2.2 x 10~ 5 


0.00015 


0.0011 


0.0053 




LOFARII-3yr 


marginalized 


0.0086 


0.0033 


0.00053 


0.0021 


0.018 (0.0050) 


0.056 


3600 


+ galaxies + Planck 


unmarginalizcd 


0.00034 


0.00030 


2.1 x 10~ 5 


0.00014 


0.0011 


0.0052 




SKA-3yr 


marginalized 
unmarginalizcd 


0.029 
0.00020 


0.011 
0.00023 


0.0016 
1.3 x 10~ 5 


0.013 
0.00024 


0.044 (0.026) 
0.0018 


0.12 
0.011 


310 


SKA-3yr 


marginalized 


0.012 


0.0036 


0.00068 


0.0027 


0.036 (0.011) 


0.12 


760 


+ Planck 


unmarginalizcd 


0.00020 


0.00023 


1.3 x 10~ 5 


0.00022 


0.0017 


0.0097 




SKA-3yr 


marginalized 


0.017 


0.0029 


0.00096 


0.0022 


0.018 (0.0031) 


0.054 


6000 


+ galaxies 


unmarginalizcd 


0.00018 


0.00019 


1.1 x 10~ 5 


0.00012 


0.00083 


0.0039 




SKA-3yr 


marginalized 


0.0074 


0.0027 


0.00046 


0.0017 


0.014 (0.0020) 


0.046 


11000 


+ galaxies + Planck 


unmarginalized 


0.00018 


0.00019 


1.1 x 10~ 5 


0.00012 


0.00082 


0.0039 




galaxy lensing 


marginalized 
unmarginalizcd 


0.045 
0.00034 


0.0052 
0.00030 


0.0024 
2.2 X 10~ 5 


0.0036 
0.00015 


0.030 (0.010) 
0.0011 


0.10 
0.0053 


970 



the Friedmann and linear growth equations to calculate its 
consequences for a(z) and D B (z). 

The Fisher matrix for the uncertainties in the esti- 
mated values of Di = D B (zi) is highly non-diagonal, show- 
ing that these parameters cannot all be estimated indepen- 
dently. To understand what properties of the growth func- 
tion can actually be constrained, it is useful to decompose 
D g (z) into modes that diagonalize the Fisher matrix, as in 
equation (|20p . Figures [7] and [5] show the well constrained 
modes in such decompositions, together with forecasts of 
their signal-to-noise ratios for surveys corresponding to our 
LOFARII-3yr and SKA-3yr cases. We choose the z; to be 
evenly spaced in ln(l + z) because we find this to give the 
best results with the fewest points, but we have tried sev- 
eral other spacings and the results are, as expected, nearly 
independent of this choice. 

In the LOFARII-3yr case (figure [7]), estimates for two 
modes are forecasted to have signal-to-noise above 2. The 
best measured mode can be interpreted as the normalization 
of the power spectrum at z ~ 2, while the higher order 
modes measure the rate , acceleration and higher derivatives 
of structure growth over the redshift range, z ~ 0.5 to 4. 

In the SKA-3yr case (figure [HJ, five modes should have 
amplitude estimates with signal-to-noise greater than 1. It 
should thus be possible to map the growth function for 
z £ 0.5. In standard ACDM D g (z) is characterized by only 
two parameters, n ma ttcr and Q,a, while four or five could be 
measured in this way, providing a test of the consistency of 
this simple model and an opportunity to discover something 



tion reconstruction, there is no longer any sample variance, 
since this part of the "noise" is now the signal. The pro- 
cedure for doing such reconstructions is described in sec- 
tion \4A\ but we will not investigate it further in this paper 
since it does not probe the background cosmology directly. 



7.4 Cosmological parameters 

Within any particular model for dark energy, the standard 
cosmological parameters can be constrained directly, as de- 
scribed in section 14.21 We here investigate constraints on 
a set of six parameters {In .4, n„, V, £2a, w , w a }- A is the 
dimensionless normalization of the primordial power spec- 
trum and n s is its spec tral index. F is the p ower spec- 
trum shape parameter (lEfstathiou et al.l Il992h given by 
T ~ D. m hexp(-2n b ) (|Peacock fc Doddsl Il994h . The first 
three parameters only affect the shape and amplitude of the 
matter power spectrum. The other three determine the ex- 
pansion history a(z) and the linear growth function D g (z). 
In this parameterization the lensing results are independent 
of H (although the CMB constraints are not) and, as a re- 
sult, many of the parameter degeneracies which can cause 
numerical instabilities when inverting the Fisher matrix are 
avoided. We assume the Universe to be flat, Q m + Qa = 1. 
This parameter set conforms to that used by the Dark En- 
ergy Task Force (DETFjf) in their investigation of the rel- 
ative merits of different observational probes of dark en- 
ergy. The dark energy equation of state parameter is here 
assumed to evolve as w(a) = p/p = w + (1 — a)w a , where 



The complete 3-dimensional density field could also, in 
principle, be reconstructed. Unlike the case of growth func- 



http: / / www. nsf.gov / mps / ast / dctf . j sp 
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Table 3. As Tabic [2] except that the sky coverage of the 21 cm lensing surveys is taken to be 25% and that of the galaxy lensing survey 
to be 48%. 







In A 


n s 


r 






Wo 


w a 


FOM 




assumed values 


-9.6 (erg = 0.75) 


1 


0.19 




0.3 


-1.0 







LOFARII-3yr 


marginalized 


0.052 


0.0068 


0.0028 




0.0049 


0.041 (0.013) 


0.13 


570 


+ galaxies 


unmarginalizcd 


0.00049 
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marginalized 


0.029 


0.0046 


0.0016 




0.0036 


0.029 (0.0048) 


0.091 


2300 


+ galaxies 


unmarginalizcd 


0.00032 


0.00032 


2.0 x 10- 
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0.00019 


0.0013 


0.0064 




SKA-3yr 


marginalized 


0.0098 


0.0035 


0.00061 




0.0022 


0.020 (0.0028) 


0.064 


5600 


+ galaxis + Planck 


unmarginalizcd 


0.00032 


0.00032 


2.0 x 10" 


5 


0.00018 


0.0013 


0.0062 




galaxies 


marginalized 


0.065 


0.0075 


0.0034 




0.0053 


0.043 (0.014) 


0.15 


470 




unmarginalizcd 


0.00049 


0.00044 


3.2 x 10" 
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0.00021 


0.0016 


0.0077 





a = (1 + ■ This is an arbitrary, but useful model that 
has become popular in the literature. 

Table[2]shows forecasts of the uncertainties in estimates 
of these parameters both unmarginalized and marginalized 
over all other parameters in the set. The difference between 
the marginalized and unmarginalized values indicates the 
importance of parameter degeneracies. For ease of compari- 
son, these forecasts (unrealistically) assume full sky coverage 
for the 21 cm and galaxy lensing surveys; the quoted uncer- 
tainties scale as f~f}J 2 for lower sky coverage accept when 
CMB constraints are included. The first four parameters are 
constrained even in the LOFARII-3yr case, but the same is 
not true for the parameters of the dark energy equation of 
state. This reflects the fact that, in the model we are assum- 
ing, dark energy is insignificant for the growth of structure 
at z it 1, while, as demonstrated in section [7^31 21 cm lensing 
is most sensitive to structure growth above z ~ 1. 

The DETF introduced a figure of merit (FOM) in order 
to compare the power of different observational techniques 
to constrain dark energy. The figure of merit is inversely 
proportional to the area of the error ellipse in the w -w a 
plane after marginalizing over the other parameters Q These 
values are given in the last column of Table [2] 

The parameter degeneracies in 21 cm lensing can be 
significantly reduced by incorporating information at low 
redshift from galaxy lensing surveys. The noise in power 
spectrum estimates from such surveys can be written as 
N K (£) = a^/n g where n g is the angular number density of 
background galaxies and er e is the root-mean-square intrinsic 
ellipticity of those galaxies. This neglects all systematic er- 
rors as well as photometric redshift uncertainties. The latter 
can be important for tomographic measurements. Following 
standard assumptions, we model the redshift distribution 
of usable galaxies as r\{z) oc z 2 e~ ( - z ^ z °' 1 , where z is set 
by the desired median redshift, and we adopt a t — 0.25. 
The proposed satellite SNAFjf] is expected to achieve a us- 
able galaxy density of n g ~ 100 arcmin" 2 with a median 
redshift z ~ 1.23, but it would only survey ~ 2% of the 

7 The figure of merit is normalized in different ways by different 
authors, and even within the DETF report itself, which defines it 
as the inverse of the area within the two sigma error ellipse but 
usually quotes it as 47r times this, as we do here. The absolute 
value has no particular significance so we will stick with what 
seems to be the de facto DETF convention. 

8 snap.lbl.gov 



sky at this depth. The DUNrfl mission proposes to survey 
almost the whole extragalactic sky to a usable density of 
n g ~ 35 arcmin" 2 with a median redshift of z ~ 0.9. Several 
planned ground-based surveys - LSStQ PanSTARRlfE 
VfST^B 

will cover comparable areas to DUNE at a similar 
depth. Since the larger sky coverage more than outweighs 
the reduced depth when constraining cosmological parame- 
ters, we adopt the "DUNE" parameter set when comparing 
the power of 21 cm, galaxy, and combined lensing surveys. In 
order to use tomographic information, we divide the galax- 
ies into 10 redshift bins each containing the same number of 
galaxies. 

We also include the predicted constraints from CMB ob- 
servations with the Planck Surveyor Satellite. This is done 
by adding the Fisher matrix for the l ensing survey s to th e 
Planck Fisher matrix as calculated in iRassat et al.l (|200cf ) . 
The CMB alone puts almost no constraint on w and w a , 
but the parameter degeneracies in the lensing surveys are 
greatly reduced by including CMB information. Our lens- 
ing Fisher matrix c o mput er code was tested against that of 
lAmara fc Refregierl (|2007l ) in order to ensure the accuracy 
of the uncertainties we forecast for parameter estimates. 

Table f2] compares the uncertainties forecast for a 21 cm 
lensing survey alone with those forecast for the galaxy lens- 
ing survey alone, and for combinations of surveys. A full-sky 
21 cm survey with LOFARII-3yr parameters would not con- 
strain these parameters as well as a full-sky "DUNE" sur- 
vey. The full-sky SKA-3yr survey would produce constraints 
that are similar to the galaxy survey on this parameter set. 
As noted above, all the action occurs at relatively low red- 
shift for this particular dark energy model. When the two 
types of surveys are combined they give much tighter con- 
straints than either survey alone. Even the seemingly insen- 
sitive LOFARII-like observations can improve the dark en- 
ergy constraints (a factor of 1.3 in the FOM) when combined 
with the galaxy survey. Combining a survey with SKA-like 
sensitivity with a galaxy lensing survey improves the FOM 
by a factor of ~ 6 and when CMB data is included the 
improvement is more than on order of magnitude. 

In practice, 21 cm lensing surveys are not likely to be 
full-sky, whereas at least half the sky could be covered by a 

9 www.dune-mission.net 

10 www.lsst.org 

11 pan-stars.ifa.hawaii.edu 

12 www.vista.ac.uk 
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satellite-based galaxy lensing survey (regions at low Galactic 
latitude are unusable). In Table[3]we combine 21 cm lensing 
surveys covering 25% of the sky with a galaxy lensing survey 
covering half the sky and the Planck CMB constraints. Even 
with this reduced size the combination of the two surveys 
results in a drastic improvement in the dark energy FOM 
and in other parameter constraints when compared with any 
of the surveys by itself. 

It is important to recognize that the constraints on dark 
energy found here for 21 cm surveys alone are in large part 
a consequence of the particular parameterization adopted 
for the dark energy equation of state. It could be that the 
energy density of dark energy is still significant above z ~ 1 
or that the rate of structure formation is otherwise signifi- 
cantly modified at such redshifts. This is the case in many 
more physicall y motivated models, for example, early dar k 
energy mode ls (jRatra fc Peeble s 1988; Caldwell et al.|[T998h . 
4-D gravity |Dvali et al.ll200Cf ). coupled scalar field theories 
(|Amendolall2000l ) and som e possible modificat ions of the un- 
derlying theory of gravity |Carroll et al.|[20o3 . for example) . 
In some of these models the growth factor could differ from 
the standard one by a factor of 2 at z = 2. Galaxy lensing 
alone would not able to distinguish between many of these 
theories, because of the poor leverage it offers on the high 
redshift Universe. As shown in section 17.31 21 cm lensing 
could probe the energy density evolution back to z ~ 7. 



8 CONCLUSIONS 

We have shown that a very large amount of cosmological in- 
formation could be extracted from the gravitational lensing 
of pregalactic 21 cm radiation. The planned low-frequency 
radio telescopes, an upgraded LOFAR or SKA, will already 
have sufficient sensitivity to do this successfully as long as 
reionization occurs sufficiently late, foreground contamina- 
tion can be accurately removed, and enough sky can be ob- 
served. Under these conditions, the matter power spectrum 
should be well measured in two and three dimensions, and 
the growth of structure will be probed from z = 0.5 to z = 7. 
Such observations would constrain many cosmological quan- 
tities, in particular the form of the primordial power spec- 
trum and presence of early dark energy, to unprecedented 
accuracy. 

Combining 21 cm lensing surveys with surveys of lens- 
ing of foreground galaxies would dramatically improve con- 
straints on structure growth at z ^ 1 where the strongest 
effects occur in the most popular parameterizations of dark 
energy. Increasing the collecting area of the telescopes by a 
factor of ~ 2 would also dramatically increase their sensi- 
tivity to lensing, as would extending their frequency range 
to lower values, if foregrounds permit. 
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